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1.  INTRODUCTION 


Experimental  testing  and  gasdynamic  modeling  of  the  ram  acceleration  technique  for  in-bore  projectile 
propulsion  are  being  investigated  at  the  U.S.  Army  Research  Laboratory  (ARL)  under  the  Hybrid  In-Bore 
RAM  (HIRAM)  propulsion  program  (Kruczynski  1994a,  1994b;  Nusca  1994a,  1994b;  Liberatore  1994). 
This  research  program  seeks  to  provide  a  highly  efficient  method  of  achieving  hypervelocity  km/s) 
projectile  gun-launch  for  use  in  high-speed  impact  testing  applications.  The  ARL  ram  accelerator  system 
uses  a  120-mm  (bore  diameter)  tube  that  is  modeled  after  the  38-mm  system  at  the  University  of 
Washington  (Hertzberg,  Brackner,  and  Bogdanoff  1988),  where  the  technology  was  first  demonstrated. 
Ram  acceleration  technology  has  also  been  successfixlly  demonstrated  at  the  Institute  of  St  Louis  (ISL) 
in  France  (Giraud,  Legendre,  and  Simon  1991;  Smeets  et  al.  1994). 

Numerical  solutions  of  the  Navier-Stokes  equations  have  been  obtained  at  the  ARL  via  computational 

fluid  dynamics  (CFD)  for  nonreacting  and  reacting,  two-  and  three-dimensional  (2-D  and  3-D)  flows. 

These  codes  are  being  used  to  investigate  the  complex  gasdynamic  physics  of  ram  accelerator  projectile 

propulsion.  A  variety  of  CFD  techniques  have  been  brought  to  bear  on  this  problem,  including  models 

for  chemically  frozen  (nonreacting)  gas,  finite-rate  global  and  multiple  step  chemical  kinetics,  and 

equilibrium  chemical  processes.  Accurate  numerical  simulation  of  hydrocarbon-based  reacting  flow  creates 

♦ 

a  very  great  demand  on  computational  resources  since  the  number  of  intermediate  species  and  the  number 
of  kinetic  steps  for  typical  hydrocarbon  fuels  are  prohibitively  large.  Global  reaction  mechanisms  based 
on  up  to  three  steps  have  therefore  been  investigated  for  use  in  preliminary  design  studies.  Using 
supercomputers,  viscous  and  chemically  reacting  gasdynamic  simulations  can  be  used  within  hours  to 
assess  the  influence  of  projectile  velocity,  tube  fill  pressure  and  mixture  composition,  and  projectile 
geometry  on  species  consumption,  mbe  wall  pressure  and  projectile  thrust  These  smdies  are  being  used 
to  seek  optimum  performance  for  the  ARL  ram  accelerator  with  minimal  gun  firings. 

The  ARL  ram  accelerator  facility,  projectile,  operation,  and  experimental  results  have  been  reviewed 
by  Kruczynski  and  are  not  repeated  in  this  report  (see  Kruczynski  1994a). 

2.  REACTING  FLOW  MODEL 

CFD  flow  simulations  for  the  ARL  ram  accelerator  are  performed  using  the  Rockwell  Science  Caiter 
Unified  Solution  Algorithm  Real  Gas  (USA-RG)  code  (Chakravarthy  et  al.  1985;  Palaniswamy  and 
Chakravarthy  1989;  Palaniswamy,  Ota,  and  Chakravarthy  1991).  This  CFD  code  solves  the  full,  3-D, 
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unsteady  Reynolds-Averaged  Navier-Stokes  (RANS)  equations  including  equations  for  chemical  kinetics 
(finite-iate  and  equilibrium).  These  partial  differential  equations  are  cast  in  conservation  form  and 
converted  to  algebraic  equations  using  an  upwind  finite-volume  formulation.  Solution  takes  place  on  a 
mesh  of  nodes  distributed  in  a  zonal  fashion  around  the  projectile  and  throughout  the  flow  field  such  that 
sharp  geometric  comers  and  other  details  are  accurately  represented.  The  conservation  law  form  of  the 
equations  assures  that  the  end  states  of  regions  of  discontinuity  (shocks,  detonations,  deflagrations)  are 
physically  correct  even  when  smeared  over  a  few  computational  cells.  The  Total  Variation  Diminishing 
(TVD)  technique  is  employed  to  discretize  inertia  terms  of  the  conservation  equations,  while  the  viscous 
terms  are  evaluated  using  an  unbiased  stencil.  Flux  computations  across  cell  boundaries  are  based  on 
Roe’s  (1981)  scheme  for  hyperbolic  equations.  Spatial  accuracy  of  third-order  can  be  maintained  in 
regions  of  the  flow  field  with  continuous  variation,  while  slope  limiting,  used  near  large  flow  gradirats, 
reduces  the  accuracy  locally  to  avoid  spurious  oscillations. 

Previous  investigators  have  used  CFD  to  simulate  ram  accelerator  flow  fields.  In  each  case,  CFD  was 
used  to  solve  the  Navier-Stokes  equations  for  the  chemically  reacting  flow  field  in  the  accelerator  tube 
aroimd  the  projectile.  IXfferences  in  CFD  analysis  arise  from  choices  of  numerical  solution/discretization 
algorithm,  representatirai  of  the  turbulence  terms,  equation  of  state,  and  chemical  kinetics  scteme. 
Soetrisno  (Soetrisno,  hnlay,  and  Roberts  1992)  used  the  LU-SGS  implicit  finite-volume  method 
(flux-vector  splitting  and  symmetric  TVD)  with  an  j^)proximate  diagonal  chemical  source  term  Jacobian, 
algebraic  turbulence  model,  perfect-gas  equation  of  state,  and  a  quasi-global,  9-species/12-reaction, 
methane-oxygen  finite-rate  kinetics  model.  Yungster  (Yungster  and  Rabinowitz  1994)  used  die  LU-SSOR 
implicit  finite-volume  method  (implicit  factorization  and  symmetric  TVD)  with  a  full  chemical  source  term 
Jacobian,  algebraic  turbulence  model,  perfect-gas  equation  of  state,  and  a  20-species/52-reaction,  finite-rate 
kinetics  model.  Li  (Li  et  al.  1992)  used  the  Flux-Corrected  Transport  (FCT)  explicit  finite-volume  method 
with  a  noninertial  source  term  to  account  for  projectile  acceleration,  inviscid  flow  model,  perfect-gas 
equation  of  state,  and  a  simplified  two-step  parametric  hydrogen-air  kinetics  model  (induction  step 
followed  by  energy  release  step).  Hosangadi  (Hosangadi  et  aL  1993)  used  the  Roe/TVD  upwind  implicit 
finite-volume  method  (Riemann  scheme  modified  for  multiphase  flows)  with  a  fully-coupled  chemical 
source  term  Jacobian,  two-equation  and  large-eddy  turbulence  models,  virial  equation  of  state,  and  a 
quasi-global,  9-species/12-reaction  finite-rate  kinetics  model. 

The  RANS  equations  for  2-D/axisyrmnetric  reacting  flow  (N  species  mixture)  are  written  in  the 
following  conservation  form  (Palaniswamy,  Ota,  and  Chakravarthy  1991): 
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The  shear  stress  terms  are  given  by 
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The  species  viscosity  and  themial  conductivity  are  referenced  to  ^nd  using  Sutherland’s  law, 
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where  Tp  and  S  vary  with  species  (StuU  and  Prophet  1971).  The  mixture  viscosity  and  thermal 
conductivity  are  determined  using  Wilke’s  law  (Wilke  1950)  denoting  f  as  p  or  k, 
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Pick’s  law  is  used  to  relate  the  mixture  diffusivity  to  the  mixture  viscosity  through  die  Schmidt  number 
Sc  =  p„/(pD).  The  specific  heat  and  enthalpy  of  each  species  (per  mass)  are  given  by  the  following 
fourth-order  polynomial  curve  fits  in  terms  of  temperature,  T,  (Drummond,  Rogers,  and  Hussaini  1987) 
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The  mixture  enthalpy,  total  energy  per  unit  volume,  and  ratio  of  specific  heats  are  given  by 


■>  -  E?.l  "i J\.«  *  E?. ,  e  =  P  *  E?.l  CT 


Y  =  1  + 


pm  _  1 


(8) 


2.1  Equation  of  State.  There  are  a  number  of  formulas  that  have  been  proposed  to  describe  the 
thermodynamic  behavior  for  states  where  the  deviations  from  a  perfect  gas  are  distinct  (Graschka  and 
Wecken  1971).  In  most  cases,  the  improvement  starts  from  the  perfect  gas  concept  as  the  limiting  case 
for  small  pressures  and  large  specific  volumes,  adapting  the  formulas  according  to  the  observed  behavior 
by  means  of  stiitable  corrections. 
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Some  of  the  best  known  examples  are  given  below. 


(p  +  _L)  (V  -  b)  =  nR„T  (van  der  Waals)  (10) 
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p(V  -  b)  =  nR„T; 
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Equations  10  and  11  are  the  simplest  ways  of  taking  into  account  die  real  gas  effects  using  molecular 
attraction  (a)  and  covolume  (b).  Equation  11  can  be  derived  from  the  Virial  equation  (Equation  12) 
assniming  veiy  high  temperatures  and  gas  densities  low  enough  that  the  mean  volume  available  per 
molecule  is  large  compared  to  the  volume  of  the  hard-sj^ere  core  of  the  molecule  (Reif  1965). 
Equation  12  is  very  general  and  in  principle  exactty  valid.  The  convergence  of  the  series,  however,  is 
guaranteed  only  for  moderately  dense  gases.  For  high  densities  in  chemically  reacting  flows,  the  formula 
is  practically  of  no  use  (Gruschka  and  Wecken  1971)  unless  the  series  is  truncated  and  certain  assumptions 
are  made  for  B  and  C.  Expressions  for  B  and  C,  the  second  and  third  Virial  coefficients,  are  given  by 
(Smith  and  Van  Ness  1987) 
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For  a  binary  mixture,  these  expressions  reduce  to 

B  =  OjBjj  +  2OJO2B22  +  02®22’  ^  ~  ^^1^2^112  3^1^2^122  ®2^22 

Second  and  third  pure  coefficients  (Bjj,  B22.  C797),  as  well  as  cross-coefficients  (Bj2.  C112,  Cj22). 

are  determined  from  experimental  data  that  is,  in  general,  available  only  for  very  limited  applications 
(Smith  and  Van  Ness  1987).  Although  Bjj^  and  are  determinable  in  principle  from  intermolecular 
potential  functions,  these  are  rarely  known  accurately  except  for  the  simplest  of  molecules.  The  virial 
coefficients  for  pure  species  can  be  obtained  experimentally  from  volumetric  data,  but  those  for  species 
in  solution  are  rarely  available.  It  is  thus  usually  necessary  to  postulate  mixing  rules  relating  the 
coefficients  of  the  solution  to  the  pure  species.  The  main  disadvantage  of  the  virial  equation  is  its 
inapplicability  to  high  densities  (Smith  and  Missen  1982).  Hosangadi  (Hosangadi  et  al.  1993)  used  the 
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Virial  equation  for  ram  accelerator  flow  fields  and  showed  some  effect  in  the  location  and  strength  of 
shock  wave  reflections.  In  these  computations  (for  a  multicomponent  hydrocarbon  mixture),  the 
cross-coefficients  for  B  and  all  coefficients  for  C  were  neglected;  thus,  C  =  0  and  B  =  SjOjBj.  For  the 
computations  described  in  the  present  work,  a  covolume  equation  of  state  is  used  (Equation  11)  where  the 
covolume,  b,  for  each  species  is  obtained  from  Hirschfelder  (Hirshfelder,  Curtis,  and  Bird  1954). 


2.2  Turbulence  Models.  The  system  of  Navier-Stokes  equations  is  valid  for  the  laminar  flow  of  a 
viscous,  Newtonian  fluid.  In  reality,  the  flow  will  remain  laminar  up  to  a  certain  critical  value  of  the 
Reynolds  number,  pU  L/p,  where  U  and  L  are  representative  values  of  the  velocity  and  length  scales  for 
the  considered  flow  system.  Above  this  critical  value  the  flow  becomes  turbulent  and  is  characterized  by 
the  appearance  of  fluctuations  in  all  variables  (velocity,  pressure,  density,  temperature,  etc.)  around  mean 
values.  These  fluctuations  are  statistical  in  nature  and  cannot  be  described  in  a  detenninistic  way. 
However,  they  can  be  computed  numerically  using  approaches  such  as  large  eddy  simulations  (LES)  and 
Reynolds-Averaged  Navier-Stokes  (RANS)  simulations. 


LES  uses  a  fixed  spacial  resolution,  with  the  effects  of  unresolved  turbulent  eddies  modeled  using 
gradient  transport  mediods  such  as  eddy  viscosity.  RANS  models  all  turbulent  fluctuations  theoretically 
or  empirically,  not  just  the  ones  smaller  than  the  computational  grid  spacing.  The  RANS  equations  are 
derived  by  decomposing  the  dependent  variables  in  the  conservation  equations  into  time  mean  (obtained 
over  an  appropriate  time  interval)  and  fluctuating  components  and  then  time  averaging  the  entire  equatiorL 
This  time  averaging  is  defined  in  such  a  way  as  to  remove  the  influence  of  the  turbulent  fluctuations  while 
not  destroying  the  time  dependence  associated  with  other  time-dependent  phenomena  with  time  scales 
distinct  from  those  of  turbulence.  After  substitution  of  these  variables  into  flie  mass,  momentum,  and 
energy  conservation  equations  and  canceling  terms  that  are  identically  zero,  two  types  of  stress  terms  are 
produced — laminar-like  stress  gradients  for  the  mean  motion  and  apparent  stress  gradients  due  to  the 
transport  of  momentum  by  turbulent  fluctuations.  These  apparent  stresses  are  commonly  called  the 
Reynolds  stresses.  For  first-order  closure  models,  the  Reynolds  stresses  are  expressed  in  terms  of  an  eddy 
viscosity,  resulting  from  the  Boissineq  assumption  (Anderson,  Tarmehill,  and  Plethcher  1984).  The 
Baldwin-Lomax  model  (Baldwin  and  Lomax  1978)  describes  the  eddy  viscosity  using  y^^  as  the  value 
where  r  attains  the  value  intermittency,  F,  and  constants  a  =  0.3,  p  =  1.6,  A  =  26. 


Pt  =  00168pFy„^r^;  F  = 


VI 


1  +  5.5 


ay 


^ymaxy 


;  r=y(l I 


(14) 
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Where  y^.  =  (yp^)(T^gjj/p)“^^.  This  model  has  been  used  extensively  for  thin  attached  shear  layers  at 
moderate  Mach  numbers  with  acceptable  results.  However,  when  separation  of  the  boundary  layer  is 
approached,  poor  predictions  result  To  improve  the  predictive  capability  of  separated  flows  using  RANS 
equations,  a  new  turbulence  model  has  been  developed  (Goldberg  1986).  The  model  is  based  on 
experimental  observations  of  detached  flows.  It  prescribes  turbulence  kinetic  energy  and  dissipation 
analytically  within  back  flows.  A  formula  for  the  eddy  viscosity  distribution  within  back  flows  is  derived 
and  used  for  the  RANS  equations  when  the  calculations  are  done  inside  separation  bubbles.  Outside  of 
them,  the  Baldwin-Lomax  turbulence  model  (Baldwin  and  Lomax  1978)  supplies  the  values  of  eddy 
viscosity.  While  the  Baldwin-Lomax  turbulence  model  is  used  to  detect  flow  separation  and  to  initiate 
ai^lication  of  the  backflow  model,  the  backflow  model  can  relocate  the  separation  point. 

2.3  Chemical  Reactions.  Hydrocarbon  reactions  are  commonly  used  for  ram  accelerator  testing  at 
ARL  (Kruczynski  1994a).  For  example. 


3CH4  +  2O2  +  ION2 . 


(15) 


This  reaction  is  one  of  a  general  class  of  hydrocarbon  reactions  investigated  by  Westbrook  and  Dryer 
(Westbrook  and  Dryer  1981) 


CaHbOe  +  1 


a  + 


b 

2 


+  c 


bo 


aCO  +  JiHoO , 
2  ^ 


(16) 


The  luel  equivalence  ratio,  <|>,  is  a  measure  of  the  amount  of  fuel  available  for  combustion.  When  (|>  =  1, 
the  proportions  of  fuel-to-oxygen  are  stoichiometric.  For  <|>  values  smaller  than  unity,  fuel  should  be 
completely  consumed.  For  the  reaction  stated  above,  (|>  =  3.0,  which  is  considered  outside  tire  range  of 
well  imderstood  CH4/O2  chemical  kinetics,  especially  for  P  >  10  atm  (Anderson  and  Kotlar  1991). 
Including  all  intermediate  species  and  reaction  steps  for  Equation  IS,  within  tiie  framework  of  a  CFD 
calculaticm  could  be  defeated  by  uncertainties  in  the  thermodynamic,  transport,  and  chemical-kinetic 
properties  of  these  species  especially  at  high  pressures  where  these  data  have  not  been  measured,  hi 
addition,  the  computational  cost  of  a  given  reaction  mechanism  depends  primarily  on  the  number  of 


8 


chemical  species  included,  rather  than  the  number  of  reactions  (Westbrook  and  Dryer  1981).  Such 
extensive  kinetic  mechanisms  (12—19  species)  have  been  included  in  CFD  analyses  using  atmosi^eric 
pressure  kinetics  data  (Yungster  and  Rabinowitz  1994).  Given  the  above  concerns,  global  reaction 
mechanisms  for  hydrocarbon  fuel  and  oxidizer  (ne^ecting  intermediate  steps)  have  been  proposed 
(Westbrook  and  Dryer  1981,  1984).  Using  one-step  chemistry. 


CH4  +  2O2  =>  CO2  +  2H2O 


(17) 


or  assuming  three-step  chemistry. 


CH4  +  2O2  =>  CO  +  2H2  +  1.5O2; 


CO  +  .5O2  =>  CO2;  2H2  +  O2  =>  2H2O. 


(18) 


The  reaction  rate  is  defined  using  the  Law  of  Mass  Action  and  an  Arrhenius  expression  for  C,  the  specific 
reaction  rate  constant: 


(0  =  Cnf.i  ol"  =  AT  “exp 


RuT 

V  “  / 


a  b  c  d  e  f 
^CH4^02^C0%2^C02%20 


(19) 


Here,  AT“  is  the  collision  frequency  and  the  exponential  term  is  the  Boltzmann  factor.  The  exponents 
a  through  f  are  given  in  Table  1. 


Table  1.  Reaction  Rate  Equation  Data 


Reaction 

Ea 

A 

a 

b 

B 

B 

e 

f 

Eqn.  17 

30.0 

8.3  X  10^ 

-0.3 

1.3 

0.0 

0.0 

1.0 

2.0 

Eqn.  18,  (1) 

30.0 

8.3  X  10^ 

1.0 

2.0 

1.0 

2.0 

0.0 

0.0 

Eqn.  18,  (2) 

44.7 

3.0  X  l(f 

0.0 

0.5 

1.0 

0.0 

1.0 

0.0 

Eqn.  18,  (3) 

37.6 

1.0  X  10^ 

0.0 

1.0 

0.0 

2.0 

0.0 

2.0 
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3.  COMPUTATIONAL  GRID 


The  governing  equations,  in  finite-volume  fonn,  are  solved  on  a  finite  mesh  of  nodes  distributed  in 
a  zonal  fashion  around  the  projectile  and  throughout  the  flow  field.  Using  this  technique,  shaip  geometric 
comers  and  contour  details  arc  accurately  represented.  Figure  1  shows  the  projectile  surface  and  tube 
surface  for  a  3-D  mesh  around  the  projectile.  The  numbers  of  axial,  radial,  and  azimuthal  grid  nodes  are 
200,  25,  and  66,  respectively.  Separate  zones  were  placed  ahead  of  the  projectile  nosetip,  over  the 
forebody,  between  the  fins  on  the  afterbody,  and  in  the  wake  regioa  The  mesh  for  an  axisymmetric 
computation  is  simpler  since  the  fins  are  removed;  therefore,  requiring  only  one  azimuthal  plane.  Separate 
zones  are  generated  as  discussed  previously.  Typically,  the  numbers  of  axial  and  radial  grid  nodes  are 
165  and  100,  respectively. 


Figure  1.  3-D  computational  grid  over  surface  of  projectile  with  four  fins  tone  hidden  from  view)  and 
over  surface  of  accelerator  tube. 

4.  COMPUTATIONAL  RESULTS 

Computational  simulations  for  the  ARL  120-nun  ram  accelerator  were  performed  to  validate  the  CFD 
code  and  investigate  various  gasdynamic  and  chemical  kinetic  effects.  Results  are  presented  in 
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sections  4.1  through  4.3.  Some  limited  3-D  numerical  simulations  are  presented  in  section  4.1.  CFD  code 
validation  results  for  axisymmetric  flow  are  presented  in  section  4.2  where  measured  and  computed  tube 
wall  pressures  are  compared.  Time-accurate  solutions  are  used  to  predict  projectile  performance  and  are 
compared  to  measured  data  in  section  4.3. 

4.1  Steady  Three-Dimensional.  Figure  2  shows  the  computed  temperature  contours  on  the  projectile 
surface  and  accelerator  tube  surface  for  frozen  (i.e.,  nonreacting)  flow.  The  Mach  number,  tube  charge 
pressure,  and  Reynolds  number  are  3.5, 51  atm  and  4.3  million,  respectively.  In  these  gray-scale  figures, 
the  lighter  colors  correspond  to  hot  temperatures.  The  forebody-afterbody  junction  is  illustrated  by  a 
vertical  black  line  (0.5L).  Three  fins  are  visible — ^top  and  bottom  fins  shown  in  planform  and  the  third 
fin  shown  on  edge  (fin  edges  illustrated  by  two  horizontal  black  lines).  Note  the  hot  gas  temperature  on 
the  fin  edge  at  the  tip  as  well  as  on  the  tube  surface  where  fin  shocks  intersect  (whitish  "X"  patterns  in 
Figure  2).  Gas  temperatures  on  the  forebody  (cone)  surface  are  uniform  except  where  the  nose  shock 
intersects  the  projectile  after  reflecting  from  the  tube  waU.  Several  expansion  regions  and  shock 
reflections  can  be  observed  along  the  fin  planform  (Figure  2,  top  and  bottom);  a  dark  (cool)  expansion 
region  over  the  fin  bevel,  a  bright  (hot)  shock  impingement  region  at  mid-planform,  and  a  succession  of 
weakening  (gray  to  black  color)  shock  reflections  along  the  fin.  A  dark  (cool)  flow  expansion  is  also 
observed  on  the  projectile  afterbody  between  fins.  The  fin-fin  shock  interaction  seen  in  the  CFD 
simulation  emulates  experimental  observations  (see  Kruczynski  1994a  where  an  experimental  flow 
visualization  via  smear  photograph  shows  combustion  emanating  from  the  fin  leading  edges  and  engulfing 
the  projectile  afterbody). 

Figure  3  shows  pressure  distributions  on  the  tube  waU  over  the  fin  tip  and  between  fins.  High 
pressures  at  the  fin  tip  (x/L  =  0.65),  between  the  fin  and  the  tube  waU  (0.6-mm  gap),  as  weU  as  between 
the  fins  due  to  shock  reflections  and  fin-shock  intersections,  are  observed.  Also  shown  in  Figure  3  is  a 
comparison  between  pressure  computed  on  the  tube  waU  for  a  3-D  computation  (between  fins)  and  an 
axisymmetric  computation  (projectile  without  fins).  Pressure  peaks  are  observed  in  similar  locations,  but 
of  reduced  magnitude  for  the  axisymmetric  case.  These  viscous,  frozen  chemistry,  3-D  computations 
required  10  CPU  hours  (CRAY  2).  Computations  with  finite-rate  (three-step)  kinetics  required  15-20  CPU 
hours  per  case. 
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1.760  K.  Black  =  394  K. 


Figure  3.  Computed  tube  wall  pressures,  Mach  number  =  3.5.  fill  p  =  51  atm.  Axisymmetric  (narrow 
line).  3-D  between  Fms  (fliick  line),  and  3-D  over  fin  top  (broken  line). 
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Much  can  be  learned  by  neglecting  the  projectile  fins  and  investigating  certain  aspects  of  ram 
accelerator  propulsion  using  axisymmetric  simulations.  These  are  the  subject  of  the  remainder  of  this 
section  while  further  3-D  simulations  wUl  be  addressed  in  future  studies. 

4.2  Steady  Axisymmetric.  CFD  simulations  were  performed  for  conditions  corresponding  to  shots  14, 
15, 22, 23,  and  25  of  the  ARL  120-mm  ram  accelerator  (Kruczynski  1994a).  While  the  gas  mixture  used 
for  shot  14  consisted  of  nitrogen  alone  (i.e.,  inert  shot),  the  other  shots  cited  used  the  mixture  described 
by  Equation  15.  The  projectile  was  accelerated  to  approximately  1,200  m/s  in  each  shot  before  entrance 
to  the  ram  accelerator  tube.  This  tube  was  pressurized  to  750,  1,000,  1,250,  and  1,500  psi  (51,  68,  85, 
and  102  atm),  for  shots  14/15, 22, 23,  and  25,  respectively  (gas  temperature  of  298  K).  Pressure  taps  were 
installed  at  0.3,  0.6,  1.15,  2.3,  3.46,  4.06,  and  4.35  m  fi-om  the  accelerator  entrance  diaphragm  and  will 
be  hereafter  referred  to  as  stations  1-7.  The  velocity  at  each  station  was  measured  using  a  radar  technique 
(Kruczynski  1994a). 

Figures  4-8  show  the  comparison  of  measured  and  computed  tube  wall  pressure  distributions  above 
the  projectile  at  a  particular  tube  station  and  for  each  shot  (see  Nusca  1994a  for  further  results).  Figure  4 
shows  the  results  for  inert  (chemically  frozen)  shot  14  using  nitrogen  gas.  As  noted  in  the  comparison 
of  3-D  and  axisymmetric  simulations  (Figure  3),  locations  of  the  pressure  peaks  agree  with  measurements, 
however,  at  a  reduced  magnitude.  The  trend  in  pressure  peak  magnitude  with  x/L  is  similar  for 
corriputation  and  experiment.  It  should  be  noted  that  the  obturator  (not  included  in  the  computation)  is 
in  relatively  close  proximity  to  the  projectile  and  could  be  affecting  these  pressures  (shot  14  only).  The 
relatively  gradual  obturator  discard  is  a  characteristic  of  inert  shots  not  observed  in  shots  through  reacting 
mixtures  (as  determined  by  monitoring  the  obturator  position  relative  to  the  projectile  base  at  shot  exit 
using  a  high-speed  camera). 

Figure  5  shows  data  for  station  4  and  shot  15,  the  lowest  pressure  shot.  The  first  substantial  pressure 
rise  (x/L  =  0.55)  is  slightly  (5%)  overpredicted.  This  j)eak  is  due  to  the  first  impact  of  the  nose  shock 
reflected  from  the  tube  wall  to  the  projectile  surface  and  back  to  the  tube  wall  again  Computations 
indicate  some  level  of  reaction  at  this  point  (Nusca  1994a).  The  second  predominant  peak  occurs  nearly 
at  the  end  of  the  projectile  (x/L  =1)  and  is  either  predicted  in  the  proper  location  or  rearward;  the 
magnitude  is,  in  some  cases,  slightly  over  or  under  predicted.  Pressures  downstream  of  the  projectile 
(x/L  >  1)  are  generally  computed  near  the  measured  magnitude,  but  without  high-frequency  oscillations. 
Several  theories  concerning  the  absence  of  pressure  oscillations  in  these  computations  are  under 
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investigation,  including  the  time-averaging  technique,  turbulence  model  employed,  and  the  assumption  of 
a  rigid  accelerator  tube. 

Figures  6  and  7  show  data  and  computations  for  shots  22  and  23,  each  representing  subsequently 
higher  tube  fill  pressures.  Results  for  station  4  (middle  of  the  tube)  are  shown  and  can  be  compared  to 
station  4  (Figure  5)  for  the  low-pressure  shot.  The  first  major  pressure  rise  (x/L  =  0.55)  is  preceded  by 
a  smaller  pressure  spike,  to  some  degree  replicated  by  the  computation.  More  pressure  spikes  are  recoided 
between  x/L  =  0.55  and  x/L  =  1  compared  to  shot  15.  As  with  shot  15,  large  pressure  rises  occur  near 
the  end  of  the  projectile  (x/L  =  1)  and  are  simulated  to  a  better  degree  for  shot  22  (Figure  6). 


Figure  6.  Tube  wall  pressures,  shot  22,  station  4.  Mach  number  =  3.44,  fill  p  =  68  atm.  Measured 
(narrow  line)  and  computed  (thick  linel. 

Figure  8  represents  shot  25  using  the  highest  tube  fill  pressure.  Comparing  this  data  (for  station  1) 
with  that  for  all  other  shots  (for  station  4)  we  find  a  trend  toward  similar  magnitude  for  the  first  pressure 
spike,  located  near  x/L  =  0.55,  and  the  second  spike,  located  near  x/L  =  1.  Both  locations  contain  reacting 
flow  and  indicate  combustion  moving  onto  the  forebody.  Data  recorded  for  subsequent  stations  (Nusca 
1994a)  for  shot  25  showed  an  "unstart"  (i.e.,  combustion  ahead  of  the  projectile,  see  Kruczynski  1994a). 


15 


25 


Figure  7.  Tube  wall  pressures,  shot  23,  station  4,  Mach  number  =  3.72,  fill  p  =  85  atm.  Measured 
(narrow  line)  and  computed  (thick  line). 


Figure  8.  Tube  wall  pressures,  shot  25.  station  1,  Mach  number  =  3.32.  fill  p  =  102  atm.  Measured 
(narrow  line)  and  computed  (thick  line). 
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The  discussion  and  conclusions  previously  presented  are  supported  by  flow  field  visualization 
accomplished  by  plotting  computed  temperature  fields  for  each  shot.  Figures  9  and  10  correspond  to 
conditions  for  shot  14  station  5  (Mach  =  3.3)  and  shot  15  station  4  (Mach  =  3.45),  respectively.  In  these 
gray-scale  figures,  the  lighter  colors  correspond  to  cool  temperatures.  Figure  9,  inert  shot,  shows  the 
projectile  nose  shock  reflecting  from  the  tube  waU,  impinging  near  the  projectile  forebody/aftetbody 
junction,  and  subsequently  reflecting  and  weakening  downstream.  High-temperature  boundary  layers  are 
observed  on  the  projectile  and  tube  walls.  Figure  10,  combustion  shot,  shows  wide  regions  of  high- 
temperature  gas  caused  by  shocks  and  combustion,  superimposed.  Note,  combustion  occurs  in  the 
projectile  and  tube  wall  boundary  layers  as  weU;  rapid  boundary  layer  thickening  caused  by  shock  wave 
impingement  is  evident.  Some  localized  combustion  occurs  in  the  forebody  boimdary  layer.  Results  for 
higher  pressure  shots  (Nusca  1994b)  show  that  extensive  combustion  between  the  projectile  forebody  and 
tube  waU  becomes  more  predominant  with  increasing  freestream  pressure.  In  addition,  increased  projectile 
acceleration  causes  the  nose  shock  to  reflect  further  rearward  on  the  tube  waU,  combustion  to  move  from 
the  second  shock  reflection  to  the  first,  and  combustion  in  the  forebody  boundary  layer  to  move  forward. 
These  effects  shift  high  pressure  to  the  forebody  where  drag  is  generated,  counteracting  thrust  generated 
by  high  pressure  on  the  afterbody. 


Figure  9.  Gray-scale  temperature  contours,  shot  14,  station  5,  inert  flow,  Mach  number  =  3.3.  fiU 
p  =  51  atm.  White  =  298  K.  Black  =  417  K. 


Figure  10.  Gray-scale  temperature  contours,  shot  15.  station  4,  Mach  number  =  3.45,  fiU  p  =  51  atm. 
White  =  298  K.  Black  =  745-2,700  K. 
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4.3  Unsteady  Axisvmmctric.  Time-accurate  CFD  simulations  are  reported  for  conditions 
corresponding  to  shot  27  of  the  ARL  120-mm  ram  accelerator  (Kruczynski  1994a).  The  projectile  was 
injected  at  approximately  1,250  m/s  (Mach  number  =  3.5)  into  the  accelerator  tube  filled  with  a  gaseous 
mixture  (see  Equation  15)  at  300  K  and  pressurized  to  750  psi  (50  atm).  The  mixture  has  a  sound  speed 
of  361  m/s,  a  Chapman- Jouget  detonation  speed  of  1,448  m/s,  and  a  precombustion  y  of  1.379  (see 
Liberatore  1994).  Pressures  on  the  accelerator  tube  wall  were  collected  at  1 1  ports  located  from  0.3  to 
9.1  m  from  the  entrance  diaphragm.  Doppler  radar  was  used  to  obtain  a  continuous  in-bore  velocity-time 
history  of  the  projectile.  The  CFD  simulation  was  started  with  the  entrance  of  the  projectile/obturator 
combination  into  the  accelerator  tube  at  t  =  20.05  ms  with  specified  velocity  (1,256.4  m/s)  and  obturator 
back  pressure  (4,000  psi).  Obturator  back  pressure  was  measured  using  a  pressure  guage  as  the  projectile 
entered  the  accelerator.  The  projectile  fins  were  neglected  so  that  axisymmetric  simulations  could  be  used. 
The  flow  field  equations  (Equation  1)  were  then  solved  in  a  time-accurate  fashion  along  with  a  force 
computation  for  the  projectile  and  obturator,  individually.  The  velocity  of  the  projectile  and  the  relative 
separation  between  the  projectile  base  and  obturator  were  updated  for  each  time  step  (0.2  ms)  using  the 
computed,  time-dependent  force  and  given  masses.  The  backpressure  on  the  obturator  due  to  the 
conventional  charge  was  assumed  to  be  ambient  after  20.45  ms.  When  the  obturator  is  greater  than  five 
projectile  lengths  downstream  of  the  projectile  base,  its  influence  on  the  projectile  flow  field  is  negligible. 
A  downstream  outflow  condition  is  then  prescribed  one  projectile  body  length  behind  the  projectile. 

Figures  1 1  and  12  show  the  computed  time  sequence  of  projectile  obturator  separation  in  terms  of  gas 
HjO  mass  fraction  and  flow  field  Mach  number,  respectively.  Water  is  a  major  product  for  the  reaction 
given  by  Equations  17  and  18,  thus  illustrating  regions  of  significant  combustion  activity.  For  Figure  11, 
white  corresponds  to  the  absence  of  H2O  in  the  flow  field,  while  dark-grey  corresponds  to  a  water  mass 
fraction  of  0.1.  For  Figure  12,  black  corresponds  to  subsonic  flow.  Flow  stagnation  on  the  obturator, 
when  in  close  proximity  to  the  projectile,  causes  both  a  normal  shock  and  a  combustion  front  to  occur  on 
the  projectile.  The  combustion  front  is  located  on  the  projectile  forebody  while  the  normal  shock  is 
located  on  the  afterbody.  How  between  the  combustion  front  and  the  normal  shock  is  slightly  supersonic. 
As  the  obturator  is  pushed  further  back  from  the  projectile  (the  last  frame  of  Figures  11  and  12 
corresponds  to  an  obturator  location  of  1.47  m  or  2.81  projectile  lengths,  behind  the  projectile  base),  the 
nornial  shock  moves  rearward  and  is  eventually  positioned  behind  the  projectile.  At  this  point,  the 
projectile  is  in  fully  supersonic  flight  Combustion  occurs  at  the  projectile-afterbody  junction  due  to  shock 
heating,  in  tte  projectile  and  tube  wall  boundary  layers  due  to  stagnation  heating,  and  downstream  of  the 
normal  shock  traveling  behind  the  projectile.  This  normal  shock,  which  also  initiates  combustion,  is  now 
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Figure  11.  Flow  field  HoO  mass  fraction  contours  (0  <  o  <  0.1).  Standard  mixture  at  51  atm.  Time 
sequence  from  20.05  ms  to  21.05  ms.  0.2-ms  interval.  Computed  projectile  velocities:  1.256. 
1,267,  1,289,  1,299,  1,305,  and  1.310  m/s.  Computed  obturator  locations:  0,  0.433,  0.961, 
1.51.  2.1,  and  2.81  m  fiom  projectile  base. 
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Figure  12.  Flow  field  Mach  number  contours  (0.09  <  Mach  number  <  3.9).  Standard  mixture  at  51  atm. 

Time  sequence  from  20.05  ms  to  21.05  ms,  0.2-nis  interval.  Computed  projectile  velocities: 
1.256,  1.267,  1,289.  1.299,  1,305,  and  1.310  m/s.  Computed  obturator  locations:  0,  0.433, 
0.961,  1.51,  2.1,  and  2.81  m  from  projectile  base. 


20 


driven  by  the  obturator  (i.e.,  a  piston).  Due  to  the  high  combustion  pressures  on  the  side  faring  the 
projectile  as  well  as  the  relieved  backpressure,  the  obturator  is  moving  slower  that  the  projectile  (this 
computational  simulation  can  be  compared  to  the  experimentally  photographed  sequence,  see  Kruczynski 
1994a). 

Figure  13  shows  the  measured  projectile  velocity-time  history.  The  projectile  velocity  was  measured 
using  a  Doppler  radar  positioned  outside  the  accelerator  and  beamed  down  the  tube  bore  using  a  sacrificial 
mirror.  Tube-mounted  transducers  were  used  to  measure  the  projectile  pressure  signature  as  it  passed  each 
station.  The  computer  simulation  shows  a  jump  in  the  velocity  relative  to  measurements  just  beyond 
station  2.  Thereafter  the  computed  projectile  velocity  was  fairly  constant  at  2.3%  above  measurements. 
The  large  projectile  thrust  computed  near  station  2  is  due  to  a  higher  than  measured  pressure  on  the 
projectile  afterbody.  It  is  thought  that  peihaps  the  obturator  is  tilted  (i.e.,  the  obturator  length  is  such  that 
the  obturator  can  spin  within  the  tube  cross  section),  thus  relieving  some  pressure  from  behind  the 
projectile.  In  the  computations,  the  obturator  is  constrained  to  move  perpendicular  to  the  projectile, 
whereas  obturator  tilting  has  been  observed  in  experiments  (Kruczynski  1994b). 


Figure  13.  computed  and  measured  projectile  velocity  vs.  time  of  flight  for  shot  27.  Mixture: 

3CH^  +20o  -H10N2  at  51  atm  and  300  K  fift  conditions. 

Figure  14  shows  a  comparison  of  the  measured  and  computed  tube  waU  pressure  distribution  at 
station  2.  Figure  14  also  shows  results  for  station  9  where  better  agreement  is  obtained.  A  possible  cause 
for  discrepancy  in  computed  pressures  are  flow  transients  introduced  by  diaphragm  rapture  and  the 
aforementioned  obturator  tilting  that  are  not  included  in  the  simulation. 
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Figure  14.  Computed  and  measured  tube  wall  pressure  distributions  for  shot  27,  station  numbers  2  and  9. 
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5.  CONCLUSIONS 


Computational  fluid  dynamics  simulations  of  the  ram  accelerator  show  that  shock  waves  originating 
from  the  projectile  nosetip  and  reflecting  from  the  accelerator  tube  wall  provide  sufficient  heating  to 
sustain  combustion  of  the  gas  mixture.  WaU  boundary  layers  on  the  projectile  forebody  and  afterbody 
sections,  as  well  as  on  the  accelerator  tube,  provide  locally  reacting  regions.  On  the  afterbody,  combustion 
is  generally  confined  to  these  wall  layers  that  thicken  with  distance.  Unsteady  simulations  show  stagnation 
of  the  combustible  gas  on  the  projectile  obmrator,  which  causes  formation  of  a  normal  shock  and  a 
combustion  wave  in  the  flow  field.  As  the  obturator  is  gasdynamicaUy  discarded,  the  normal  shock  trails 
behind  the  projectile  and  the  combustion  wave  collapses  into  shock-induced  combustion  on  the  projectile. 
This  effect  has  been  observed  computationally  as  well  as  with  experimental  flow  visualization.  With  the 
obturator  sufficiently  downstream,  the  projectile  accelerates  at  supersonic  speeds.  Projectile  velocity 
through  the  accelerator  can  be  predicted. 
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LIST  OF  SYMBOLS 


b  covolume 

Cp  specific  heat  capacity,  constant  p 

Cy  specific  heat  capacity,  constant  volume 

C  specific  reaction  rate  constant 

D  mass  diffusion  coefficient 

e  specific  total  internal  energy 

Eg  activation  energy 

F,  G  flux  vectors  (Equation  1) 

g  Gibbs  energy  per  mole 

h  molar  specific  enthalpy 

i  i-th  species 

L  total  length  of  projectile 

m  mixture  quantity 

M  molecular  weight 

n  number  of  moles 

N  total  number  of  species 

p  static  pressure 

Pr  Prandtl  number 

Q  heat  content  of  the  mixture 

R  specific  gas  constant,  (y  -  IjCp/y 
Rg  universal  gas  constant,  RM^ 

Sc  Schmidt  Number,  p^/pD 

t  time 

T  static  temperature 

u  axial  velocity 

V  radial  velocity 

U  magnitude  of  the  local  velocity  vector 

V  specific  volume  (1/p) 

W  dependent  variable  vector  (Equation  1) 
X,  y  cartesian  coordinates 

X  species  mole  fraction 
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Y  ratio  of  specific  heats,  Cp/c^ 

AHf  enthalpy  of  formation 

K  heat  transfer  coefficient 

p  molecular  viscosity 

V  stoichiometric  coefficient 

p  density 

a  species  mass  fraction 

%  shear  stress  tensor 

©  chemical  production  term  (Equation  19) 

Q  source  term  vector  (Equation  1) 
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